Scale-dependent bias from primordial non-Gaussianity in general relativity 
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In this note we examine the derivation of scale-dependent bias due to primordial non-Gaussianity 
of the local type in the context of general relativity. We justify the use of the Poisson equation 
in general relativistic perturbation theory and thus the derivation of scale-dependent bias as a test 
of primordial non-Gaussianity, using the spherical collapse model. The corollary is that the form 
of scale-dependent bias does not receive general relativistic corrections on scales larger than the 
Hubble radius. This leads to a formally divergent correlation function for biased tracers of the mass 
distribution which we discuss. 
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I. INTRODUCTION 

One of the most powerful discriminants between vari- 
ous models for the origin of structure in the Universe is 
the non-Gaussianity of the primordial fluctuations. The 
non-linear evolution of inhomogeneities on super-Hubble 
scales during or after slow-roll inflation can give rise to 
the local form of non-Gaussianity where the primordial 
curvature perturbation can be given as a local function 
of a Gaussian random field P, 0, H, 0] . This occurs, for 
example in the curvaton scenario, due to the late-decay 
of a weakly coupled scalar field after inflation or in 
the new ekpyrotic models of an accelerated collapse with 
two or more fields @, 0, 0, H - 

Local non-Gaussianity is commonly characterised by a 
single parameter, /nl, where the local primordial New- 
tonian potential on large scales in the matter-dominated 
era is given by [1 01 ] 

Ma;) = 0g(z) + /nl (<&(» - (<t>%)) , (1) 

and 4>g( x ) is a Gaussian random field 1 . 

In Newtonian gravity the Poisson equation, 



>N 



= AirGpS , 



(2) 



gives the local density contrast, S, as a function not of 
the local Newtonian potential directly, but rather as a 



1 Note that we adopt the sign convention of Komatsu & Spergel 
|10( which is widely used to characterise primordial non- 
Gaussianity; the Newtonian potential is positive for a point mass 
</>N( r ) = +Gm/r. Many theoretical calculations of primordial 
non-Gaussianity, starting with that of Maldacena [ill , use the 
opposite sign for the Newtonian potential, as for instance used 
in the influential review by Mukhanov, Feldman and Branden- 
berger ll'J . and this leads to a different sign for /nl- Although 
the choice of sign is conventional, observations are sensitive to 
the sign of / NL . 



function of the spatial divergence of the Newtonian po- 
tential. A consequence of combining equation ([T]) with 
the Poisson equation ([2]) is that the non-Gaussian part 
of the density field is not a local function of the Gaussian 
part. 

If we re-write the physical Laplacian in terms of co- 
moving coordinates in a Friedmann-Robertson- Walker 
(FEW) universe, V 2 = d 2 /a 2 , where a is the scale factor, 
we can write the Poission equation as 



-d 2 



>N 



= 4TrGa 2 pS . 



The Fourier transform of this equation gives 
fc 2 Mk) = ^ 2 H 2 S(k) 



(3) 



(4) 



where we use the Friedmann equation to give the back- 
ground density in terms of the Hubble expansion, H. It is 
the scale-dependence of this relation which is ultimately 
responsible for the diverging scale-dependent bias due to 
non-Gaussianity as k — > 0. 

It is commonly assumed that the formation of astro- 
physical objects is a local function of the dark matter 
overdensity (smoothed over some scale R): 

<5 9 (x) = /(5(x), (V<5(x)) 2 , V 2 <5(x), . . .) + S (x), (5) 

where s(x) is a non-deterministic stochastic component. 
Iff R is large enough so that 5 is small and if Fourier 
modes are decoupled, then the overdensity of astrophysi- 
cal objects traces dark matter fluctuations in the Fourier 
space in the limit of k — > 0: 



S g (k) =M(k) + «(k), 



(0) 



where b = /'(0) is the so-called bias parameter. The 
corresponding power spectrum is given by 



P gg (k)=b 2 P(k) +l 



(T) 



The stochastic white noise tail has never been observed 
and hence we will neglect it here. The important point is 
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that bias b should approach a constant as k approaches 
zero. Deviation from this implies either that the for- 
mation process is non-local or that Fourier modes are 
coupled, i.e., the dark matter density fluctuations non- 
Gaussian. 

In concrete models of formation of astrophysical ob- 
jects, it is assumed that the function / is such, that it 
puts objects in the peaks of the local density field [HI]. 
Such density peaks collapse under their own gravity and 
form virialized objects. The local number density of such 
objects in Lagrangian space is given by 



n = no(l + bhSi), 



(8) 



where stands for Lagrangian space bias, Si is the con- 
tribution from the very long wavelength modes that es- 
sentially modulate the mean density of the effective local 
cosmology. Therefore 



bh =n — 



_i dn 
dSi 



(9) 



and the more usual Eulerian-space bias is given by b — 
1 + &l (Til. [HI HH [ri| . This approach towards formation 
of astrophysical objects is the essence of the so-called 
peak-background split. If one understands how the local 
effective cosmology is modulated by the presence of the 
large-scale modes and how the formation of objects is af- 
fected by the change in the local effective cosmology, then 
one can calculate the bias parameter for that particular 
tracer. 

The bias is therefore given by the interplay between 
the small-scale collapse dynamics and the large scales, 
which modulate the local cosmology through fluctuations 
present in the very large scale modes. The purpose of 
this paper is to analyse how can this be viewed within 
the context of general relativity. 

The distinctive scale-dependence of large-scale struc- 
ture has recently been proposed as a powerful test of 
the local form of primordial non-Gaussianity 18]. For a 
biased tracer of structure, it can be shown that local non- 
Gaussianity induces a scale dependent bias b = be + Aft. 
with [H: 



where 



Ab 



a(k) 



2/nl dlogn 
a(k) 9 log erg 



2c 2 T(k)D(z) 



(10) 



i Hi 



is the Fourier space conversion factor between the present 
day density perturbations and the Newtonian poten- 
tial during matter domination on large scales, including 
the linear transfer function T(fc), which is approaches 
unity on large scales, and growth factor normalised 
to be D(z) = 1/(1 + z) in the matter domination. 
d log n/d log as describes how the number density of as- 
trophysical objects, n, is affected by the change of the 



amplitude of small-scale fluctuations as in Gaussian cos- 
mologies. Therefore, although the initial derivation relied 
on the statistics of Gaussian fields, it is now clear that 
the effect is present for all local theories for the formation 
of the astrophysical objects [2(| [2l|, [HI, H3- In particu- 
lar, for a universal form of the mass function n(M), the 
above equation simplifies to 



A6 = / NL (6 G -1)- 



i H 



c 2 k 2 T(k)D(z) 



(12) 



where be is the Gaussian bias and <5* ~ 1.68 is the linear 
over-density at collapse for the spherical collapse model 
and other symbols have their usual meaning. 

A striking feature of the Equation (Tl2"|) is that the 
scale-dependent correction diverges as the wave-vector 
k approaches zero. This is a direct consequence of the 
Newtonian form of the Poisson equation in the Fourier 
space. Since the Newtonian approximation is expected 
to break down at the scales comparable to or larger than 
the Hubble scale, it is timely to ask whether the diver- 
gence of the scale-dependent bias is real or an artefact of 
the Newtonian approximation. This is the main question 
that we address in this brief report. 

We start by analysing the Poisson equation in the con- 
text of general relativistic perturbation theory in Section 
HI! Next we review the well known exact general rela- 
tivistic solution, namely the spherical collapse solution 
f Section IIIip . By investigating the metric in the limit of 
small perturbations, we can identify the gauge in which 
the spherhical collapse calculation is performed. This 
justifies the form the Poisson equation for linear pertur- 
bations on all scales. In Section IIVI we discuss the for- 
mally divergent correlation function of a biased tracer of 
structure and show how this is absent in observations of 
real correlation functions. We conclude in Section M 



II. GENERAL RELATIVISTIC POISSON 
EQUATION 

One might naively expect corrections to the Newto- 
nian Poisson equation of order (aH/k) 2 which would 
become large close to the Hubble scale and change the 
scale-dependence of a(k) on large scales. However, the 
form of the general relativistic constraint equations is 
gauge dependent. 

In a completely homogeneous universe, the FRW met- 



( U ) ric ds 2 



dt 2 



(dn 2 



8ijdx 1 dx :l , 



has 



no gauge freedom. At the first order of perturbations, on 
the other hand, the coordinate freedom gives rise to am- 
biguities in the definition of the density perturbation in 
an inhomogeneous spacetime (see, for instance, Ref. 24]). 
In particular the freedom to make an inhomogeneous re- 
definition of the time coordinate n — ► n+Sn(x, n) leads to 
a redefinition of the density perturbation Sp — > Sp — p'5n, 
where a prime denotes derivatives with respect to confor- 
mal time, rj. 
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The most general scalar perturbation 2 of the spatially 
flat Friedman-Robertson- Walker (FRW) metric can be 
written as 0, E| 

ds 2 = a 2 {(l + 2A)dr) 2 - 2(d l B)dx i dr 1 

- [(1 - 2i/>)6ij + 2(d i d j E)]dx i dx j } . (13) 

but the metric perturbations A, B, ip and E are all, like 
the density perturbation, 5p, gauge-dependent. 

In an arbitrary gauge there is a first-order energy con- 
straint 

-m (ip' + HA)+d 2 [ip + H{E' - B)} = AirGa 2 Sp . (14) 

where TL = aH is the conformal Hubble rate, and a first- 
order momentum constraint 



di (ip' + HA) = -A-KGa 2 di{5q) 



(15) 



where di (Sq) is the 3-momentum. We can combine these 
two equations to obtain a general relativistic Poisson 
equation 

d 2 [tp + H{E' - B)] = AttGo 2 [Sp - 3H6q] . (16) 

This is exactly the same as the Newtonian Poisson 
equation ([3]) as long as we identify —<pN with the metric 
potential in the Newtonian or longitudinal gauge 



>N 



= = V + H{E' - B) , 



(17) 



and p8 with the density perturbation in the comoving- 
orthogonal gauge 



pS = Sp c — Sp — ZTLSq . 



(18) 
= 



i.e., the density perturbation in a gauge in which B 
and 8q = 0. 

The general relativistic Poisson equation (Til)]) relates 
gauge-invariant combinations of metric and matter per- 
turbations, which nonetheless have the interpretation as 
potential or density perturbations in two different gauges 

iii si. 

If the primordial Newtonian potential, (f>pf, has a lo- 
cal non-Gaussianity as given in Equation (JTJ) then the 
comoving-orthogonal density contrast, 5, will indeed have 
a non-local non-Gaussianity. 

Note however, that in some gauges there are general 
relativistic corrections to the Poisson equation. For ex- 
ample, in the longitudinal gauge the shear potential, 
E 1 — B, vanishes and on large scales the regular solution, 
— 4>n — Ai — ipi = const, gives a relativistic correction 
to the Poisson equation 



Thus on super- Hubble scales (\d 2 (j) N \ < W 2 |<^aH) the 
density contrast in the longitudinal gauge, Si — Spij p, 
is proportional to the Newtonian potential 

5t(x) ~ -2(f> N (x) 

= -2cj> G {x) - 2/ NL (0 2 (x) - (^D) . (20) 

This density contrast, Si, does have a local non- 
Gaussianity in the large scale limit. 

Both the density perturbation in the comoving- 
orthogonal and longitudinal gauge can be given as gauge- 
invariant combinations of gauge-dependent variables, and 
on sub-Hubble scales (\d 2 ipN\ ^ W 2 |?W|) we find from 
Eqs. (TiT)|) and (flTp)) that Spi ~ Sp c and these two possi- 
ble definitions of the density perturbation coincide. The 
question then arises as to which is the relevant density 
perturbation from which to calculate the bias of galaxies 
and large-scale structure. 



III. SPHERICAL COLLAPSE 

The usual Newtonian calculation of the abundance of 
collapsed objects such as galaxies or quasars assumes 
that the collapse is a function of the local density field 
smoothed on some scale R. But which is the appropri- 
ate density variable as we approach the Hubble scale? 
Although initial conditions are set at early times where 
the linear perturbation theory is valid, collapse requires 
a non-linear calculation which is generally only possible 
in general relativity if we make some assumption of sim- 
plifying symmetry. The most simple example is that of 
spherical collapse. In the absence of any pressure per- 
turbation (as for pressureless matter, with or without a 
cosmological constant) the model allows us to compute 
the non-linear collapse of a top-hat overdensity described 
by a closed FRW model embedded within a flat FRW ex- 
terior. 

The Friedman equation for the unperturbed outer uni- 
verse is given by: 



H z = [ - 

a 



SirGpo 



(21) 



while the inner region has a perturbed scale factor, a, 
and small positive curvature, K, and hence 



H 



~ \ 2 
a 



8irGpo 
3a 3 



K 

^2' 



(22) 



m 2 cj) N 



d 2 (p N = ^Ga 2 S Pi 



(19) 



Without the loss of generality, we can set 8nGpo /3 = 1 . 
The Friedman equation for the closed inner universe can 
then be expressed in a parametric form 



2 We neglect vector and tensor perturbations which are decoupled 
from scalar perturbations at first order. 



1 . 

a = ( 1 — cos ( 



t = 



1 



2^3/2 



(23) 
(24) 



4 



Expanding a and t to fifth order in 9, one can show that 
in the limit of 9 -c 1, 



a= — 



2/3 



(25) 



We see that to leading order the scale factors match as a 
approaches zero, if we have a = (3</2) 2 / 3 as the solution 
for the exterior flat universe. 

The linear over-densities are thus given by 



a 3 a 20 



-Ka, (26) 
5 



where Sa — a — a. We have thus recovered the expected 
result that linear perturbations grow with a in the matter 
era and that this spherically symmetrical system has a 
one parameter solution. When the inner perturbation 
collapses, 9* — 2n, t* = t:K~ 3 / 2 and the linear density 
at that moment is 



<5* 



2/3 



1.68, 



(27) 



which is the standard result. 

But to decide which linear density perturbation this 
is, we need to identify the coordinate choice implicit in 
the spherical collapse model. Let us look at the same 
process from the metric perturbation perspective. Inside 
the overdense region the perturbed, non-linear metric is 
given by a closed FRW metric 



ds 2 



dt 2 



h 2 



dr 2 



1 - Kr 2 



r 2 dQ? 



(28) 



where K > 0. For Kr 2 <C 1, this can be recast in Carte- 
sian form 

ds 2 = dt 2 - a 2 [5ij + Kr 2 (d l r)(d J r)] dx l dx j . (29) 

Comparing this with the linearly perturbed FRW metric 
of equation (|13[) we can identify the scalar metric pertur- 
bations: 



J, 



Kr 



4 

E : 



A = 


(30) 


B = 


(31) 


Sa 


(32) 


a 


Kr 4 


(33) 
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In addition, the spherical collapse picture is implicitly set 
up using comoving spatial coordinates, the expan- 

sion or collapse is determined solely by the local scale fac- 
tor, a, and there are no peculiar velocities. Since B = 
we are working in a comoving-orthogonal system with 
vanishing 3-momentum, Sq = 0. 

Since A = 0, we are also working in a synchronous 
coordinate system [which was implicitly assumed in the 



use of the same time derivative in Eqs. (|2Tj) and (|22|) ]. In a 
comoving-orthogonal coordinate system the conservation 
of momentum requires 



6P C 
P + P 



(34) 



For pressureless matter momentum conservation thus re- 
quires A = in a comoving-orthogonal gauge, i.e., the 
gauge is also synchronous [27] , so that the proper cosmic 
time in the overdense interior and unperturbed exterior 
are the same. 

Thus we see that the initial data for the spherical col- 
lapse problem which is set up at early times in the matter 
era where we may use relativistic perturbation theory is 
set up in a comoving-orthogonal gauge. In particular the 
relevant density perturbation which is used to determine 
whether a given region will collapse by the present time 
is the comoving-orthogonal density perturbation defined 
in Equation (fl8|) which is related to the Newtonian po- 
tential by the general relativistic Poisson equation (Ti"7o|) . 

Note that the general relativistic Poisson equation (JTSJ) 
shows that for the growing mode solution 4>n — const, we 



have Sp c 



as the physical spatial gradients decay 



in an expanding universe. Thus the density contrast, 
5 = Sp c /p oc a grows on all scales in the matter era, 
as in Newtonian gravity and in the linearised spherical 
collapse solution, Equation ([26]) . 

Since the Poisson equation, ((T5|) . relates the New- 
tonian potential (fTTJ) to the density perturbation in 
the comoving-orthogonal gauge (TTg)) . and the spheri- 
cal collapse model uses the density perturbation in the 
comoving-orthogonal gauge to determine the collapse of 
overdensities, we conclude that the Newtonian form of 
the scale-dependent bias does not receive any corrections 
arising from general relativity on the super- Hubble scales 
in the spherical collapse model. 



A. Spherical collapse in longitudinal gauge 

It is of course always possible to consider spherical col- 
lapse in different coordinates. If we redefine the confor- 
mal time and radial coordinate 

1 - ^K(r 2 - r 2 t )\ m , 

r = [l + ^I<V 2 i+K(r 2 -r 2 )y e , (35) 

then the line element (|29|) can be written (to first order 
in K) in a longitudinal gauge 

ds 2 = a 2 [(1 - 2(t> N )drij - (1 + 2(j> N ) (dr} + rjdti 2 )] , 

(36) 

where ai = a(r)e) and the Newtonian potential for ri < r$ 
is 



(37) 
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Note that we need to specify the size of the overdensity, 
7"o, so that the coordinate time coincides with the time 
in the unperturbed background, r > r$. 

Note that the density contrast in the longitudinal 
gauge, given by Eq. (TT9"j) . is no longer spatially homo- 
geneous 



St = -Ken 
5 



1 

1 + r 



' H 



(38) 



where th = T~C 1 is the comoving Hubble length. The lin- 
ear central density contrast at the collapse time is given 
by 



x i 3 f 3lr 

fc.|r=0-g( Y 



2/3 



(39) 



In contrast to the result in the comoving-orthogonal 
gauge, Eq. (f!?T|) , the criterion for collapse in the longi- 
tudinal gauge is not simply a function of the local den- 
sity, but it also depends upon the scale of the overdensity 
for tq ~ rn- Thus, although the GR correction to the 
Poisson equation in the longitudinal gauge, Eq. (|19D . im- 
plies that the density contrast on super-Hubble scales 
(ro 3> rn) becomes a local function of the Newtonian 
potential, the collapse condition, Eq. (f3"9"| . becomes non- 
local and this is the origin of the scale-dependent bias as 
seen in the longitudinal gauge. 

More generally, dropping the condition of spherical 
symmetry, we can combine the two different expressions 
of the GR Poisson equations (HHJ) and (fTO)) to obtain 
a manifestly non-local expression for the density con- 
trast in the longitudinal gauge in terms of the comoving- 
orthogonal density contrast 



St = (i-m 2 d- 



IV. CORRELATION FUNCTION 



(40) 



The necessary consequence of the fact that there are no 
general relativistic corrections to the form of the scale de- 
pendent bias on scales larger than the horizon is that the 
correlation function formally diverges as we show next. 
The correlation function of objects is given by the Fourier 
transform of the power spectrum 



£(r) 



1 



(2tt) 3 



d 3 kP(k)e 



ik-r 



(41) 



In the homogeneous and isotropic universe, this can be 
integrated to 



« r > = 5? 



p(kfj^Sllk 2 dk. (42) 
kr 



The dark matter power spectrum Pdm(fc) — k™ 3 for 
small fc, where n s is the spectral index of primordial fluc- 
tuations. Therefore, for P(k) = (b + Ab) 2 P dm (k), the 



P(k) is proportional to k n '~ i (for objects with b ^ 1) as k 
approaches zero and hence the integrand is proportional 
to k ns ~ 2 in this limit. The result is that £(r) diverges for 
every r. 

How does one then relate the measurements of the 
large scale structure to the observations in non-Gaussian 
universes? The hint is given by the fact, that root 
mean square of fluctuation in the field of tracers, i.e. 
a = y£(0) is infinite. However, when we measure the 
fluctuations in a given survey, we measure the mean den- 
sity from the same survey, effectively forcing the fluctua- 
tions to average to zero. Therefore, the measured density 
fluctuations S are given by 



<5(x) = <5(x) - 5, 



where 



/ 



d 3 xlf(x)(5(x) 



(43) 



(44) 



with IT(x) being the normalised survey window. 

This means that the measured correlation function is 
given by 



/ 



£{r) = / d 3 xIT(x)(5(x)(5(x + r) 



d 3 xVF(x)(5(x)(5(x + r) - S 2 . (45) 



The expectation value for £(r), averaged over all pos- 
sible realisations of the underlying density fields is thus 
given by 



£(r) ) = / d 3 xM/(x) (<f(x)*(x + r)) - (S 2 ) 

= t(r)-(P). (46) 



The second term of this equation is the variance of <5 
over an ensemble over possible realisations of the underly- 
ing field. As the universe is ergodic, we can replace many 
realisations of the underlying field with many copies of 
the window function over a single realisation and hence 
the second term can be written as 

(6 2 ) = ((S*W) 2 (x)) 

= 7i / d 3 kW k (k) 2 P(k) = cr 2 w , (47) 

where * denotes convolution and Wk is the Fourier trans- 
form of the window function. The last term is therefore 
nothing more than the variance of the field smoothed 
over the window function of the survey. We can there- 
fore write 



f(r))=£(r) 



(48) 
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This is an expected result. We simply cannot measure 
variance coming from scales larger than the survey. 

Let us take a concrete example. For the spherical top- 
hat window function of size R, the integral (|T7|) becomes 

a 2 (R) = At f d 3 kP(k)k 2 T 2 (kR) 7 (49) 
27T J 

with T(x) = 3 (sin x — x cos x) /x 3 . In this case, the Equa- 
tion (j4"8)) can be rewritten as 

i{r) = i / p(fc) {^hv 1 ~ T2(kR) ) k2dk - (50) 

This is now a well behaved transform. In the limit of 
k — y 0, the term in brackets becomes k 2 (R 2 /lO — r 2 /6) 
and hence the integrand is proportional to fc™ 3 and thus 
convergent. 

Note that this cancellation occurs for any finite window 
function W. If a window function has a largest typical 
linear scale R, then for k <C Wk ^ 1 as the exp(ifci?) 
is varying slowly across the window. One these scales 
s'm(kr)/kr — W 2 (kR) is therefore quadratic in k to the 
leading order. 

Finally, there are two other worries associated with a 
diverging field. First, one might worry that the are for- 
mal problems with perturbing a homogeneous universe 
with a solution that has infra-red divergence. However, 
the underlying physical field, the dark-matter density, is 
a completely regular and does not diverge on the large 
scales. Only some non-linear transformation of this reg- 
ular field, namely a biased tracer of this field, or a field 
of peak positions (after smoothing on some scale) has a 
large-scale divergence. Therefore, this should not be a 
problem. A more serious issue is that the field of tracers, 
although it can be arbitrarily biased with respect to the 
underlying dark matter density field, must still remain 
positive, i.e. 6(x) > —1 since the local density of astro- 
physical objects cannot be negative. To estimate when 
this becomes important, we demand that the dimension- 
less power per unit log k, A 2 = P(k)k 3 /An 2 is much less 
than unity on large scales. For biased tracer, this reduces 
a requirement that /nl(& — 1) <C 10 4 , which is always true 
for observationally interesting values of /nl and b. 

V. CONCLUSIONS 

In this paper we have discussed the derivation of scale- 
dependent bias from primordial non-Gaussianity in the 
context of general relativity. Our results are somewhat 
unexciting - we find no corrections arising from gen- 
eral relativistic effects on scales larger than the Hubble 
length. We have justified this by noting that the spherical 
collapse and thus peak-background split formalisms that 
are performed in the Newtonian context can be carried 



over into general relativity if one works in terms of the 
density perturbation in the comoving-orthogonal gauge. 
This justifies the use of the Newtonian form of the Pois- 
son equation ([3]) on all scales (in contrast to what one 
might naively think [28|). 

This means that the power spectrum of biased tracers 
of the density field diverges as the scales of interest be- 
come larger and larger and leads to the infinite variance 
in the local density of the biased objects. The formal 
correlation function is divergent on all scales, but any 
observed correlation function will be finite, because vari- 
ance on scales larger than the size of the survey will not 
be observed. 

The reason that general relativity reduces to appar- 
ently Newtonian equations is because we are only con- 
sidering scalar perturbations which obey energy and mo- 
mentum constraint equations relating the metric to the 
matter perturbations. Considering first-order perturba- 
tions about an FRW background, the constraint equa- 
tions enforce a Poisson-type relation between the local 
density perturbation and the metric perturbations on a 
given constant-time hypersurface, coinciding with a given 
time in the background spacetime. Coordinate freedom 
in general relativity allows different choices of time slicing 
in the perturbed spacetime, but it is the density pertur- 
bation in the comoving-orthogonal gauge that follows the 
same behaviour as the Newtonian density perturbation 
on all scales. 

Assuming spherical symmetry we are able to follow 
the non-linear collapse in a comoving-orthogonal coordi- 
nate system. Spherical symmetry eliminates the vector 
and tensor parts of the metric perturbation which might 
have introduced non-Newtonian behaviour. For example, 
there are no gravitational waves if we impose spherical 
symmetry. Deviations from spherical symmetry could 
thus introduce additional general relativistic terms, but 
the scalar, vector and tensor perturbations are decou- 
pled from scalar density perturbations at first-order, so 
such effects would only be expected to arise at second- 
or higher-order [29] . This might ultimately be important 
when we wish to constrain non-linearity parameters, /nl, 
of order unity but it suggests that general relativistic cor- 
rections to the Newtonian results currently being used are 
negligible for current data. 
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